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ABSTRACT 


During the annual Gray Whale migration from the 

۲٣٠٠٦ Sons to Baja California, the mammals travel in coastal 
waters, thereby presenting an opportunity for the study of 
their sound spectral and source levels and vocabulary.  How- 
ever, such measurements are distorted by surface and bottom 
reverberation, Using the theory of rough surface scattering, 
knowledge of the bottom impedance, and correlation techniques, 
it is possible to decompose the shallow water reverberation 
into the contributions from different paths. From this, the 
range, depth and the deverberated spectral source levels of 
the sounds of the mammal can be determined by use of only one 
hydrophone rather than the ee three or four. The 
theory, deverberation programming, and experimental results 
are presented for a model of the whale's pulsed radiation in 


a laboratory model coastal environment. 
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I, INTRODUCTION 


Once yearly, the Gray Whale, Eschrichtius glaucus, 
migrates southward from the Aleutians and passes very close 
to the California coast in shallow water. During this 
migration, the sounds in the water close to the whales can 
easily be recorded; but, they may not be the true sounds 
produced by the whale. The whale normally produces an 
intermittent, short duration signal which in shallow water 
is received at the hydrophone via direct, surface reflected 
and bottom reflected paths. Since the path lengths are 
different, the signals arrive at the hydrophone at different 
times and they interfere. To obtain the true sounds produced 
by the whales, this interference, which is called reverbera- 
tion, must be eliminated. 

The purpose of this thesis is to model in the laboratory 
the whale sounds in the shallow water and to develop a 
technique to determine the range, eliminate the surface and 
bottom reverberation and calculate the spectral source levels 


Asa function of time, 
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In a reverberant environment, the original signal can 
only be realized if the reverberation is removed. The 
method used to remove the reverberation, which is here 


t 


called '"deverberation," assumes that the whale is a point 
source, the geometrical spreading is spherical, the water is 
isovelocity and the water depth is constant. 

Before the deverberation technique can be applied, the 
direct, surface reflected and bottom reflected path dis- 
tances must be known. Normally, this information is obtained 
by knowing the source position and calculating the path 
distances from the known geometry. Determining the hori- 
zontal source position is difficult and requires at least 
two directional or three omnidirectional hydrophones, accu- 
rately fixed with respect to each other at all times. To 
determine the depth requires at least one additional hydro- 
phone at a different depth, Generally, three or four 
hydrophones are deployed [Refs. 1 and 2]. With each addi- 
tional hydrophone, the complexity of the system is increased 
Since each hydrophone requires its own amplifying and record- 
ing system. In shallow water, however, taking advantage of 
the surface and bottom reflections, one can use a single 
hydrophone and gather all the information necessary for the 
application of the deverberation technique. 

Consider the direct, surface scattered, and bottom 


scattered sounds received at only one hydrophone which 





is deployed in shallow water as depicted in figure 1. It 
will be shown that when the differences between the arrival 
times for the three paths are known the three path distances 
can be calculated. Using the surface specularly scattered 
path, Ro» it is seen that 

Rg = PNE Re 

5.557. с = 


2 
Re = С^ е! < 


where 


EC Es [R?- (HOY 


EC Eo Ls en oy) 


Substituting and rearranging terms gives 
D z 
ечи. 2 


u 5 
(2) Re sy ہہت‎ 44D) 


Using e the time difference between the direct path 
arrival and the surface reflected path arrival and C, the 


mean speed of sound, gives 
u 
C To E Rs + Re - R 


and therefore 
K: 
ect. = (R*+ 44D)” -R 
Similarly using the bottom specularly scattered path 
t | u 
Rez Rg + Rg 
2. 2 
Ве“ +۶(ہ2۰)‎ ٥ 


К" = (2-0)? + 8"? 
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where 


٢ (2-5) EMI E 5 
$ * (22-H-D) | RF (H-0) | 


CEA) 2 Е 5 
5 * Qi-u-5 ۳۰۲۰ 


Substituting and rearranging 


m ے‎ (=. ру 2 2 E 
a) 


mi- 


„ _ 2-H) 2 2; 4D - ZH-2D) 
(5) Ке“ (оз -н-5) IR + 4(2 + | 


Using t,, the time difference between the direct path 


arrival and the bottom reflected path arrival gives 
1 M 
c Ta = Кв + Re - Қ 
and, therefore, 
2 
(6) ct, = |R*4 4(z^sup-zW-zD) -R 


Solving equations (3) and (6) simultaneously for R, the 
range of the source from the receiver and D, the depth of 


the source, produces 


Mee totes 42" - чен) - T. (c Te) 
(7) О = 


4} нтр + Т. (2 - Н) | 


m 2. 
SR = —THD= (ста ) 
erh 
The equation for the range is left as a function of the water 
depth to facilitate its being programmed. Now that the range 


and depth are known, the other two path distances can be 


calculated. From equations (1) and 0۰۰۰۰۷٤۹989 80 
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Dach distance is 


J- 
(9) Re = (R* * 4u0)* 


and from equations (4) and (5), the bottom reflected path 


distance is 
n 
Q0) Rg » |R*4 4(z? sup - zu- гр)? 


Therefore, when the time arrival differences for the differ- 
ent paths are known, equations (8), (9), and (10) can be 
used to determine the first three path distances. The paths 
for multiple reflections can be calculated directly from the 
known geometry, assuming specular scatter. 

For a transient signal, determination of the differential 
arrival times To and тр сап be realized by the use of an 
autocorrelation technique. The correlation function can be 


defined as [Ref, 3] 


(11) R (t) * EI Iyw-ullvasr) -u]! 


where v is the time average, u is the mean, and E is the 
expected value of the received signal. This process is 


programmed using a digital summation 


n- 
(12) Κίτ) - a A [m0 - ul uysal 


with n representing the total number of samples in the record. 


Performing the autocorrelation on the reverberant signal at 


the one hydrophone gives peaks at delay times corresponding 
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to zero delay time, and the arrival delay times from the 
reflected signals, The peaks are realized only when the 
direct signal is delayed enough to correlate with the 
reflected signals. The computer program called AUTOPEAK 
performs the autocorrelation and then applies an envelope 
detection to determine the delay times for the peaks. These 
delay times are then used in equations (7), (8), (9), and 
(10) to determine the range, depth, and reflected path dis- 
tances. Thereby, all the geometrical information necessary 
for the deverberation technique has been obtained. 

Computer programs have been developed for deverberation 
in either the frequency or time domain. 

The computer program designed to eliminate the reverbera- 
tions in the frequency domain is called DEVERB. For the time 
being, assume only one frequency component, whose amplitude 
and frequency are functions of time. For a signal with time- 
varying frequency and amplitude, we can describe the pressure 
at the receiver, due only to the direct path signal as [Ref. 
4] 

(13) Р, (4) = TA 
Then, taking into account spherical divergence, the spatial 
phase shift, and specular scattering from a Gaussian rough 
surface, the pressure at the hydrophone due to the surface 


scattered signal can be written as 


-3 3leldr-(Rs- RARA) + Tr 
(14) RU) =P R-t) = De ac ау ] 


LIT, 
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and the pressure due to the once bottom reflected signal is 


¡[o (Ni - (Re-R)R (1) +3 ] 
(15) АЕ = Py lt - Tp) > τὰ ce 
12 Te 


R® and e 8/4 are the frequency-dependent pressure amplitude 
reflection coefficients, and Y and T. are the phase shifts 
due to the bottom and surface reflections, respectively. 
The surface reflection coefficient depends on the roughness 


of the surface. The exponent for specular scattering is 


τος. 5] 


16) 9: EL cos Os 


where 9 is the R.M,S. wave height, A is the wavelength of 
the sound signal and ©. is the angle of incidence. Equations 
(14) and (15) are derived from previously received direct 
path pressures, corrected for path differences and phase 
shifts and represent the pressures due to the single reflected 
paths, The coherent sum of equations (13), (14) and (15) is 
the pressure sensed by the receiver when these three compo- 
nents are present, 

The signal processing is done digitally. Therefore, the 
continuous time dependence is replaced, سے تم‎ it occurs 
in the previous equations, by digital block numbers, indi- 
cated by the subscript index "k."" Each block contains 
enough data samples of the incoming signal to give the desired 
spectral frequency resolution during a block duration which 
is small compared to the total duration of the time-varying 


Signal. The frequency change of the signal within a block 
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duration is assumed to be much smaller than the frequency 


resolution of the digital processing. The block duration is 


BENUMBER OF SAMPLES IN THE BLOCK 


Е SAMPLING FREQUENCY 


W i ime the block ends is related to the continuous time t 


by 


My t = TK: K=1, 2, 3 


Veti 


The index "i" is used for the spectral frequency compo- 
nent of the complex wave being analyzed. 
The equation for the source pressure at unit distance 


using the frequency deverberation correction is 


j (Ski RR) a í oe نے‎ 0× 
A 


9 : 
is wi T ARg- R)RUATY 
Fire i IDE E. 2i | : Mn) 


¿lem Cg - R)R,- 5 
کک‎ e mI PORC uy 


(18) Юк (4) e 


- 3 (“¬ -R)RL-T- YX 
[Re Die ide Зк 


where І(К-М), 1(К-М) and 1(K-L) are unity factors with values 
1(K-N) 


1 K»N 


0 otherwise 


1(К-М) = 1 K > M 


0 otherwise 
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1(K-L) 


II 
ES 


K > L 


0 otherwise 


The pressure amplitude of the 1 ΕΠ frequency component in 
block K of the receiver reverberant signal is represented 

by ki and its phase represented by PRO Ihe deverberated 
pressure amplitude is represented by De i and its phase by 
Sri The first term on the right hand side of the equation 
represents the received signal, The second term represents 
the correction due to a single specular scatter from the 
surface; the third represents a single bottom reflection 
correction, and the fourth represents the correction for a 
path which includes one surface and one bottom reflection. 
The equation can be expanded to include other multiple reflec- 
tions. 


The block indices are determined by 





Ш E 
- 

(20) Ns k - 5 
"e 

(21) L =k - ss 
τ 


The output of the frequency deverberation program is a 
series of spectra of the comsecutive blocks and its Fourier 
transform is a time plot of the deverberated signal. 

The above procedure takes the signal from the time 


domain (the time series after A/D conversion) by Fourier 
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transform to the frequency domain, where the known frequency 
dependent reflection coefficients are easily applied, and 
then back to the time domain to verify the effectiveness of 
the process. 

When the reflection coefficients can be assumed to be 
frequency-independent, a simple point-by-point deverberation 
procedure can be applied in the time domain. The applicable 


temporal deverberation equation is 
>R ye Ds etc 
y. 2 O سے‎ πως. — + 2 
22) Dk = Cys Se 2 2 Di- 8 Z, Dw + бе ¿REA 


d sample. The 


ος represents the pressure amplitude for the KE 
other terms are similar to those in equation (18). For low 
roughness surfaces g < l and the use of e78/2 over the 
appropriate frequencies will be a good approximation which 

is essentially independent of frequency. One advantage of 
the temporal deverberation technique is the relative freedom 
from restrictions of block size; the block size is determined 


only by the desired frequency resolution and the rate of 


change of frequency of the transient source. 


79 
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In order to model the Gray Whale's environment in 
Monterey Bay, a three meter cube "anechoic" water-filled 
tank was used. An artifical bottom made of hard rubber with 
an experimentally determined pc of approximately 2.4 x 10° 
mks rayls was positioned one meter below the water surface. 
This type of bottom was chosen for its specific acoustic 
impedance since the bottom at the listening area in Monterey 
Bay has a oc of approximately 3 x 109 mks rayls. The depth 
of the bottom and relative placement of the source and 
receiver were determined in order to obtain realistic modeled 
delay times between the reflected signals and the direct 
Signal. A 1.8 cm diameter spherical source was used because 
of its small size and its ability to transmit a signal with 
minimum distortion; but it also limited the minimum frequency 
to about 10 kHz, An FM up-sweep of varying widths and sweep 
rates was used to model one of the sounds produced by the 
Gray Whale. 

The equipment was connected as in figure 2 with the mas- 
ter unit pulse generator being used to synchronize the samp- 
ling frequency oscillator which determined the start and stop 
frequencies, sweep rate and pulse width. The pulse repeti- 
tion rate was determined by the master unit pulse generator. 
The signal was then amplified and sent to the source. The 
reverberent signal was received by an LC-10 hydrophone, and 


then amplified again and sent to the analog-to-digital 
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converter, The A/D converter was triggered through the 

delay generator via the unit pulse generator and the sampling 
frequency oscillator. A sampling frequency of 320 kHz was 
used, The delay generator delays the trigger by the time 
required for the signal to be transmitted through the water 
and then this delayed pulse triggers the unit pulse generator. 
The unit pulse generator was used to determine the total on- 
time of the sampling frequency oscillator. The oscillator 
determined the samples per second taken by the converter 
during the oscillator's on-time. This complex equipment set- 
up was designed to allow the A/D converter to sample only 
during the time when the direct and reflected signals were 
present, thereby reducing the amount of computer time and 
storage required to process the data. Reflections from 
other surfaces in the tank were eliminated wherever possible 
in the sampling time by varying the pulse repetition rate 

of the master pulse generator or by varying the source and 
receiver placement. After the analog signal was changed by 
the converter to digital form, it was stored on cassette 
memory and then analyzed using the programs AUTOPEAK and 


BEVERB. 
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IV. DATA PROCESSING AND RESULTS 


The autocorrelation plots of reverberation for two 
different sweep widths are seen in figure 3 with the scale 
factor for the delay time equal to 3.125 u seconds. The 
plot of the 90 kHz sweep width shows a steeper slope of the 
envelope, thus a more clearly defined peak than the one with 
only a 10 kHz sweep width, This indicates, as expected, that 
as the difference between the upper and lower frequencies 
decreases, the correlation peaks become harder to determine. 
In the limit of only one frequency being present, there would 
be no peaks in the envelope, Figure 4 shows the range error 
(computer determination compared to direct measurement) versus 
the ratio of the upper frequency to the lower one. The graph 
indicates that for a ratio above 1.2 to l the frequency sweep 
of the signal is sufficient to get accurate time differences 
for the reflections and thereby to determine the source range 
and depth from the autocorrelation processes. 

The frequency deverberation program is designed to give 
a true spectrum of the signal by eliminating frequency depen- 
dent reverberations. Since the signal is time-varying, a 
small block of time is desirable to keep the change in the 
signal to a minimum during the block. The limiting factor 
to the minimum block size is the desired frequency resolu- 
tion. For a spectrum to accurately represent an instant of 


time rather than being a spectrum averaged over a length of 
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time, the frequency change during the block time should be 
much less than the frequency resolution, In addition, where 
possible, the block duration should be equal to, or a sub- 
multiple of а - T,| im order to permit equation (18) to 

be applied, 

From figure 4, it is known that depth and range can be 
determined only when the frequency ratio of the chirp is 
greater than 1,2 to l, Figure 5 shows the reverberant and 
deverberated signals for a chirp from 46 kHz to 54 kHz with 
ratio of 1,17 to 1. The reverberant signal, top of figure 
5, was divided into blocks, equal to Tt, - Ttg, and then trans- 
formed back into the time domain which is shown at the bottom 
of the figure, Some improvement can be seen but the expected 
slowly increasing frequency at the approximately constant 
amplitude has not been realized. It is believed that the 
inadequate deverberation is due to the fact that the frequency 
sweep during a block is approximately four-tenths of the 
frequency resolution. A slower sweep rate or a larger block 
duration would have cured this problem, However, a slower 
sweep rate for the model would have decreased the accuracy of 
the range determined by the autocorrelation; and a large block 
duration was precluded by the geometry which determined 
| Tp - E The temporal deverberation technique which is 
presented next did not suffer from these limitations. 

The result of applying the temporal deverberation program 
Noa signal with a sweep width from 5 kHz to 95 kHz can be 


seen in figure 6 with the reverberant signal on the top and 
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the deverberant signal below, Since the source was resonant 
at 64 kHz, the FM sweep grows in amplitude to 64 kHz and 

then decreases. The deverberated signal shows the frequency 
and amplitude modulation cleanly until the end of the direct 
signal, The reverberation after that time is due to scatter- 
ing from the side walls of the tank. Possible reverberation 
due to the fourth and later terms on the right hand side 

of equation (22) were excluded by limiting the duration of 


sampling by the computer. 
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V, CONCLUSIONS 


The work described in this thesis demonstrates the 
feasibility of obtaining a non-reverberant spectrum of a 
transient source in a reverberant environment. The technique 
includes calcylation of the autocorrelation of the received 
signal to determine range and depth of the source and com- 
puter processing to correct for the surface and bottom reflec- 
tions, 

The autocorrelation function provides an accurate method 
for obtaining the range and depth of a source of transient 
sound in shallow water. The correlation technique can be 
performed for a chirp sound with ratio of upper to lower 
frequency of greater than 1.2 to l. At least two reflections 
are required to obtain the depth and range of the source with 
respect to the receiver, 

Frequency and time deverberation programs which use the 
position data from the autocorrelation have been developed 
to eliminate the reverberations and, thereby, to obtain 
corrected spectra or corrected time plots. Either technique 
can be used; however, because the output of the computer is 
a time series, it is natural to apply temporal deverberation. 
This becomes very simple if the surface and bottom reflection 


coefficients are independent of frequency. 
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